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Abstract 



The contribution of a 7T7r-exchange three-body force to the three-nucleon 
binding energy is calculated in terms of a ttN amplitude. The latter is based 
on a meson-theoretical model of ttN interaction developed by the Jiilich group. 
Similar to a previous study based on simple phenomenological ttN potentials 
a very small effect of the resulting three-body force is found. Possible ori- 
gins of the two-orders-of-magnitude descrepancy between the present result 
and the values obtained for the Tucson-Melbourne three-body force are in- 
vestigated. Evidence is provided that this discrepancy is most likely due to 
strikingly different off-shell properties of the ttN amplitudes underlying the 
two approaches. 
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I. INTRODUCTION 



Recently T.-Y. Saito and I. R. Afnan (SA) [|],EJ have calculated the contribution of a im- 
exchange three-body force (see Fig. [D to the three-nucleon binding energy in terms of the 
ttN amplitude using perturbation theory. Their approach determines the contributions of the 
different ttN partial waves and (via the division of the irN amplitude into a pole and nonpole 
term) allows for a consistent determination of the irNN form factor. The calculations are 
based on phenomenological separable ttN potentials H]. The total contribution of this three- 
body force (TBF) to the binding energy of the triton has been found to be very small. It is 
typically of the order of a few keV. This result falls short of calculations based on the Tucson- 
Melbourne (TM) and the Brazilian 7T7r-exchange three-nucleon potentials [[§-§]. The latter 
potentials make a contribution to the binding energy of the triton that is of the order of the 
discrepancy between experiment and calculations with realistic nucleon-nucleon potentials 
(i. e. about 1 MeV). 

The origin of this surprisingly large difference of two orders of magnitude in the contribu- 
tion of the three-body force has been the topic of two subsequent very detailed investigations 
by Saito and Afnan [[J and Murphy and Coon J7|. SA concluded that the total contribution 
of the TBF to the binding energy of the triton found in their approach is so small as a result 
of the energy dependence of the ttN amplitudes, cancellations between the contributions 
from the S- and P-wave vrA^ partial waves, and in particular, the soft nNN form factor. In- 
deed the form factors extracted from their nN interaction models correspond to monopole 
cutoff masses of around or even less than 0.4 GeV - which have to be compared to values of 
about 0.7 - 0.8 GeV suggested from other information |pHil"H and to the value of 0.8 GeV 
used in calculations with the TM potential. 

Murphy and Coon carried out a thorough comparison of the nN amplitudes that underlie 
the TM TBF and the calculations of SA f7j . They criticized that the amplitudes used by SA 
do not fulfil the low-energy theorems of the irN interaction as imposed by chiral symmetry. 
But still they attested that these amplitudes are qualitatively similar to the one used in the 
TM potential. The prime reason for the two-orders-of-magnitude discrepancy found in the 
contributions of the TBF was suspected to be likewise the soft 7rNN form factor emerging 
from the phenomenological separable ttN interactions employed by SA. 

In the present paper we want to re-investigate the origin of this large discrepancy in 
the predicted contribution of the TBF. We follow the same approach as SA. However, we 



start out from a meson-theoretical tcN model developed recently by the Jiilich group |12 
This model, besides being conceptionally much better founded than the simple separable 
potentials employed in Refs. |],|2| , has the important advantage that it does not exhibit those 
deficiencies which led to a criticism on the work by Saito and Afnan. Firstly, the Jiilich model 
is in agreement with empirical information on the ttN amplitude in the subthreshold region 



TJ|. In particular its prediction for the amplitude at the so-called Cheng-Dashen point is 



close to the empirical value. Second, and most importantly it yields irNN form factors 
which are comparable to those of the TM potential, with monopole cutoff masses of around 
0.7 GeV ||14|| . Therefore this model provides an ideal starting point for re-assessing the role 
of the 7T7r-exchange TBF in the binding of the three-nucleon system. 

The paper is structured in the following way: In Section II we review the salient features 
of the Jiilich ttN model. In particular we concentrate on those properties that are relevant 
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for the present study. The specific structure of our three-nucleon code makes it necessary 
to represent the meson-theoretical ttN amplitude in separable form. For this purpose we 
applied the so-called EST method ||15| which allows to generate separable representations 
that agree exactly (on- and half-off-shell) with the original interaction at specific predeter- 
mined energies. This method is briefly described also in Section II. Furthermore, we discuss 
the reliability of the separable representation by comparing (off-shell) amplitudes obtained 
from it to the ones of the original Jiilich model for various ttN partial waves. 

The formulation of the TBF is given in Section III together with a short outline of the 
formalism. Results for the contribution of the 7T7r-exchange TBF to the triton binding energy, 
based on the ttN amplitude of the Jiilich model are presented in Section IV. Anticipating 
our results we find again that the TBF is very small. Therefore, in Section V, we embark 
on a detailed discussion of the ttN amplitudes on which the TM TBF and our calculation 
are based. We focus specifically on the off-shell properties of these amplitudes since they 
are determined quite differently in the two approaches. Indeed, we will argue that much of 
the observed two-orders-of-magnitude discrepancy in the binding energy is due to off-shell 
effects and we will present numerical evidence for this conjecture. Finally, a summary is 
given in Section VI. 



II. THE ttN INTERACTION 

The ttN models employed in the present investigation are based on meson exchange 
and have been developed by the Jiilich Group fll2 |. They include the s-channel and u- 



channel nucleon and delta-isobar pole diagrams together with correlated tttt exchange in 
the J p = + (a) and 1~ (p) channels as shown in Fig. [| The interaction potentials are 
derived in time ordered perturbation theory and then unitarized by means of a relativistic 
(Lippmann-Schwinger type) scattering equation 

T = V + VG T . (2.1) 

The resulting models account for the scattering data in the elastic region as well as for the 



low-energy parameters fl2| . Furthermore they also satisfy chiral symmetry constraints. In 
particular, the resulting values for the so-called ttN S term (the contribution of the isoscalar 
forward scattering amplitude at the Cheng-Dashen point) - X = 66.4 (65.6) MeV for model 



1 (2) of Ref. - are in good agreement with the empirical value of £ = 60 MeV [T6 

Since in the discussion of Murphy and Coon |7j special attention was given to the sub- 
threshold behaviour of the ttN amplitudes, we want to present here the corresponding results 
for the Jiilich nN model. In the continuation of the T-matrix for the potential models to the 
subthreshold region, we follow the procedure outlined in section 4 of Ref. 0. We calculated 
the (on-shell) background isoscalar amplitude F + (u,t) (conventionally called D + ; for defi- 
nition see, e.g., Refs. |7]f|) from the subthreshold point v = 0, t = to the Cheng-Dashen 
point {y = 0, t = 2m^). The predictions of the Jiilich models 1 and 2 are shown in Fig. |3| 
in comparison to the ones of the ttN amplitudes that form the basis of the TM and Brazil 
7T7r-exchange TBF. One can see that the results of the meson-exchange models more or less 
coincide with the nN amplitude employed in the TM TBF. Furthermore, they are also in 
rather nice agreement with the empirical subthreshold amplitude given by Holder H. 
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In a recent paper C. Schiitz et al. JTj] have determined the nNN vertex functions result- 
ing from the models of Ref. 0] (cf. Ref. [|14[] and the Appendix for definitions and relevant 
formulue). From the vertex functions at the nucleon pole, the decrease in the ttNN form 
factors F nNN (p 2 ) from the pion pole p 2 = m 2 to p 2 = 0, has been extracted following the 
procedure proposed by Mizutani et al. |T7|]. This quantity is a measure for the softness of 
the nNN form factor. It has been found that the nNN form factors implied by the models 
considered in Ref. [T3j] are, in general, significantly harder than the ones used by SA in their 
study of the contribution of the 7T7r-exchange three-nucleon force to the 3N binding energy. 
In particular, the model 2' of Ref. [13"] yields a decrease in the form factor from the pion 
pole to p 2 = of slightly less then 4% - which is quite close to the value of 3% implied by 
the form factor introduced in the TM three-nucleon force and also in good agreement with a 



recent lattice QCD calculation B and other independent information flTOUTTf. Accordingly 



the major concern raised by Murphy and Coon against the work of SA does not apply for 
this model and therefore it provides an excellent starting point for re-analyzing the contri- 
butions of the 7T7r-exchange three-nucleon force to the 3N binding energy in the approach 
of Saito and Afnan. 



In this work we will also present results for the other models considered in Ref. |14 
and it is appropriate to say some words about the basic differences between these models. 
All the models are based on the same dynamical input (cf. Fig. |2|). They differ, however, 
in the (phenomenological) parametrization of the (bare) vertex form factors. Details and 
explicit formulae can be found in sect. Ill of Ref. fll2| . Both models provide a similarly good 
description of the nN scattering data. However, the differences in the parametrization of 
the vertex form factors lead to different vertex functions and in turn to different (dressed) 
nNN form factors. 

Two further models (1' and 2') have been presented in Ref. [T3J for the following reason: 
The models 1 and 2 are constructed by assuming nNN pseudovector (pv) coupling. However, 
the nNN form factor derived from lattice QCD calculations M, and the values of F v nn{0) 
from other independent information |fL0| , |lI| (and also the values given for the models used by 



SA), are based on pseudoscalar (ps) coupling. In order to allow for a meaningful comparison 



for the different couplings the authors of Ref. JI3J have constructed variants of models 1 and 



2 (labelled 1' and 2') where ps coupling is used in the nucleon s-channel pole terms. 

Values for F^nn^) f° r the various models are compiled in Table I. (Note that F^nn is 
normalized to F n jsrN(m 2 ) = 1.) For the ease of comparison we also include here the result 
for the separable nN model PJ by McLeod and Afnan || which has been used (amongst 
others) in the investigation by SA and on which most of the concerns and criticism of 
Murphy and Coon are based. It is evident that for this model the decrease of the form 
factor from the pion pole to p 2 = is much larger - almost 20%. We want to emphasize 
here, however, that it is the nNN vertex function which enters into the calculation of the 
TBF and not the form factor (cf. Section III). Therefore we also show this quantity (cf. 
Fig. |]). Note that the nNN vertex functions for models 1 and 1' are practically the same. 
Accordingly we expect these models to give the same result for the contributions of the 
7T7r-exchange TBF to the 3N binding energy, and therefore we will consider only model 1 
in the following analysis. Furthermore, the vertex function for model 2 is obviously harder 
than the one resulting from model 2' (for small and intermediate momenta) - in contrast the 
form factors for model 2' looks harder (cf. Table I). This seemingly paradoxical situation 
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has been discussed thoroughly in Ref. [14 . 



Before carrying out the actual three-nucleon calculations the meson-theoretical nN mod- 
els have to be expanded in separable form. This is necessitated by the specific structure 
of our three-nucleon code which can only deal with interaction models given in separable 
form. For this purpose we apply the so-called Ernst-Shakin-Thaler (EST) method [|TjJ which 
allows us to generate separable representations of arbitrary rank N that agree exactly (on- 
and half-off-shell) with the original reaction matrix at N specific predetermined energies. 

Let us begin with the (partial wave projected) Lippmann-Schwinger equation for the 
radial wave function 

\^ E ) = \k E ) + G (E)V\^ E ), (2.2) 



where \k E ) is the incoming wave and Go(E) the two-body Green's function. (In Eq. (|2.2| ) and 
in the following the partial wave index is suppressed for convenience.) For proper scattering 
solutions (on-shell) k E and E are related by 



E = ^m% + k% + ^ml + k%. (2.3) 

According to the EST method a rank-iV separable representation for the potential V is given 
by 

N 

V= E V^K^V , (2.4) 

where E iy (i = 1,...,N), are N freely chosen energies. The coupling strengths Ajj are 
determined by the condition 

N 

J2^(^\V\^ Ek ) = S lk . (2.5) 



It is evident from Eq. ( |2.4|) that the "form factors" of the separable potential V [18] consist 
of the objects V\4> Ei ), where \i[> Ei ) are solutions of Eq. ( [2.2|) for the potential V at the 
energies Ei. Therefore, by virtue of Eq. ( |2.5| ), the following relation holds 

V\^ Ei ) = V\^ Ei ) = T(Ei)\k Ej ) = f(Ei)\k Ei ) (2.6) 

at the N energies E iy where T is the solution of the Lippmann-Schwinger equation ( |2.1| ) for 
the separable representation V. This means that the on-shell as well as the half-off-shell 
t-matrix for both interactions V and V are exactly the same at the energies Ei. 

In the present case the interaction models V are energy-dependent and therefore a modi- 
fication of this scheme proposed by B. Pearce [O] is employed. According to it the condition 



N 

(^ El \V{E)\^ Ek ) = ^ El \V{E)\^ Ek ) = £ <^|^)|^>Ay(J5)<^l^i)IK> (2-7) 



has to be used for determining the coupling strengths Xij(E) instead of Eq. (|2.5|). As 
expected also the separable representation becomes now energy-dependent. 
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A special treatment is required for the Pu partial wave which contains the (s-channel) 
nucleon pole. It must be possible to clearly separate the contribution of this pole term from 
the total Pu amplitude. Its contribution to the three-nucleon binding energy is already 
taken into account by solving the standard bound-state Faddeev equations. Therefore, in 
order to avoid double counting, only the non-pole part of the Pu must be considered in 
the present investigation. Furthermore in the consistent approach of SA the nNN vertex 
function is extracted from this pole term and is then used for the vertices where the pions are 
emitted (absorbed) by (at) the outer nucleons. Consequently a separable representation for 
the Pu partial wave must guarantee that (a) the non-pole amplitude is reliably reproduced 
and (b) the ttNN vertex function extracted from the pole term agrees exactly with the one 
obtained for the original interaction model. This can be achieved and we summarize details 
of the construction procedure in the Appendix. 

Of course if one relies on such separable representations one has to ensure that they in- 
corporate all the relevant properties of the original interaction models. From extensive tests, 
we find that a rank-1 separable representation is sufficient for the present purpose provided 
that the expansion energy is choosen in the relevant energy domain (i. e. around the triton 
binding energy which corresponds to a center-of-mass energy of roughly 930 MeV in the ttN 
system). Thus we have selected this particular energy for the separable representation to be 
applied in the present study. Since this energy is below the elastic ttN threshold, fcg in Eq. 
( p.2|) can no longer be fixed by the on-shell condition. Following our previous work |20| we 
choose k,E in such a way that iUe fulfils Eq. (|2.3|) . For E\ = 930 MeV this implies Ue 1 ~ 138 
MeV/c. 

The quality of the separable represention is demonstrated in the Figs. [| and |] for model 
2. Fig. |5| shows the off-shell transition amplitude T a (q,q';Z) for fixed off-shell momenta 
q = q' — 130 MeV as a function of the total energy Z. Note that in case of the Pu partial 
wave only the non-pole part of the t-matrix is shown since, as explained above, this is the 
part relevant for the present study. Fig. ^| shows the transition amplitude T a (q,q'; Z) for 
fixed Z and q' as a function of the other off-shell momentum. In the latter figure we display 
only results for those partial waves that are expected to give the dominant contribution of the 
7T7r-exchange three-nucleon force to the 2>N binding energy [I], namely S31, Pu, and P 33 . We 
want to point out, however, that the agreement in the other partial waves is of similar quality. 
Likewise we want to refrain from displaying corresponding results for the model 1 here since 
the quality of its separable representation is pretty much the same. Furthermore, we do not 
show the nNN vertex functions resulting from the separable representation because - by 
construction - they are identical to the ones of the original model (cf. the Appendix). 

Note that we have also constructed rank-1 separable representations where the expansion 
energies are near the ttN threshold (E\ = 1077 MeV for Su', E>\ = 1000 MeV for S31; 
Ei = 1100 MeV for Pu, P31, P13, -P33). (In case of the S31 partial wave the amplitude 
around 1077 MeV has a peculiar energy dependence which would require, in principle, a 
rank-2 representation. In order to avoid this complication we have chosen a somewhat lower 
value for E]_.) These representations will be also employed in our investigations. They will 
serve as a term of reference for how strongly the resulting 3N binding energies depend on 
the specific choice of the separable representation. 
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III. FORMULATION OF THE THREE-BODY FORCE 



In this section we formulate the TBF using the irN amplitude and irNN vertex function 
which have been derived in the previous section. 

We follow the prescription by SA. The TBF is schematically shown in Fig. [7], i.e., (i) 
a pion is emitted from the first nucleon, (ii) the pion is scattered off the second nucleon, 
(iii) the pion is absorbed by the third nucleon. The strength function of pion emission and 
absorption is given by the ttNN vertex function which depends on the energy of the ttN 
system as explained in the previous section. Also the irN amplitude corresponding to the 
scattering of the pion on the second nucleon is energy dependent and includes only the non- 
pole contribution in the Pn channel. We introduce Jacobi variables in the irNNN system 
so that we define the ttN relative momenta and energies at the stages (i), (ii) and (iii). 

The momenta of the three nucleons before the three-body interaction in Fig. [7| are given 
by ki, k 2 and k 3 and the momenta after the interaction are k' 1; k 2 and k' 3 , respectively. We 
define the relative momenta between the third nucleon and the system of first and second 
nucleons and pion (qs), between the second nucleon and the system of first nucleon and pion 
(p 3 ) and between the first nucleon and pion (Q 3 ) at the stage (i) 

q 3 = -k 3 , (3.1) 

m N (K + ki) - (m N + m 7r )k 2 

Ps = 77; s , (3.2) 

(2m N + m w ) 

Q 3 = r — y , 3.3 

where k^ is the pion momentum. Then the center-of-mass energy in the system of first 
nucleon and pion, E 3 , is obtained using those relative momenta 

E 3 = E + mN _A_A (3.4) 
2/i 2 2/ii 

where E = —Et is the total energy of the whole system not including rest masses, and the 
reduced masses \i\ and /x 2 are defined, respectively, by the relations 



1 1 



\ii m N rriM + m^ 

and 

1 1 1 



+ ; , (3-5) 



/i 2 m N 2m N + 



(3.6) 



In the same way we define the relative momenta and energies of the system of second nucleon 
and pion, Q 3 (before scattering), Q[ (after scattering) and E 2 at the stage (ii) and Qi and 
Ei of the system of third nucleon and pion at the stage (iii), respectively. 

Using these nN relative momenta and energies the TBF, W(E), can be symbolically 
written as 

W(E) = v* N (Qi;Ei) G^NNiE) T„ N (Q[,Q' 3 ; E 2 ) G wNNN (E) v* N (Q 3 ;E 3 ). (3.7) 



7 



Here, v^ N (Qi] Ej) (i = 1, 3) is the renormalized ttNN vertex function (cf. the Appendix) 
which gives the strength of the pion emission and absorption on the nucleon. G t nnn{E) is 
the propagator of the ttNNN system and T^iQ'^Q 1 ^ E 2 ) is the non-pole part of the 11N 
scattering amplitude. Eq. (|3.7|) represents how the pion is emitted from the first nucleon, 
scattered off the second nucleon, and then absorbed by the third nucleon. Note that the ttN 
amplitude and the ttNN vertex function are determined in the same framework, i. e. they 
are obtained from the same nN interaction model. A detailed discussion of Eq. ( |3.7| ) can be 
found in Ref . @ . Since in the Jiilich model the energy is defined fully relativistic whereas a 
semi-relativistic form is employed by SA, we have to change the nNNN propagator of Ref. 
§ to be 

G„ NNN (E) = {Ei - y/Q* + m % - y/Q* + miy\ (3.8) 
with i = 1 or 3 (cf. Eq. (4.9) of Ref. @). 



IV. RESULTS 

The contribution of the TBF to the binding energy of the three-nucleon system, AE^ 3 \ 
is calculated in first order perturbation theory, i. e. 

A£ {3) = (V\W(-Er)\V), (4.1) 

where | is the triton wave function. This wave function is obtained from solving the Fad- 
deev equations for the so-called PEST potential pO| , ^I| which is a separable representation 



of the Paris NN potential derived by the EST method. All nucleon-nucleon partial 
waves with total angular momentum less than or equal two are employed in the calculation. 
The triton properties obtained by the PEST potential are comparable with those by the 
original Paris potential as shown in Refs. f22|j2|]. 

Evidently, like in the work by SA there is no consistency between the NN interaction 
which is used to generate the triton wave function and the TBF. However, we are making 
use of the Born approximation and therefore the triton wave function and the TBF are 
obtained separately, anyway. Thus, as argued already by SA , we do not expect that this 
inconsistency has an influence on the qualitative features of our results. One should also 
keep in mind that a similar inconsistency is involved in standard 3N calculations employing 
the TM TBF. 



The contributions of the TBF generated by the irN models 1, 2 and 2' of Ref. |nj 
to the 3N binding energy are listed in Table |. These results are based on the separable 
representations of the original interaction models described in the preceeding section. Note 
that the value of 930 MeV is used for the expansion energy since we expect the average 
ttN energy in the three-nucleon system to be near this value. In Table 0] also the result for 
model PJ of Ref. 0] is listed for the sake of comparison. Thus, it is easy to see that our 
results are qualitatively very similar to those of SA. Again the contibutions of the individual 
ttN partial waves to the 3N binding energy are very small (only in the order of keV) and 
there is again a strong cancellation between the attractive P\\ and P33 partial waves and 
the repulsive S31 partial wave. 
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In Ref. Murphy and Coon have presented a thorough comparison of the nN amplitudes 
used in the TM and Brazilian TBF and the one applied in the calculations by SA. Their 
main conclusion was that the two-orders-of-magnitude smaller value for the binding energy 
resulting from the 7T7r-exchange TBF based upon the separable models of SA is most likely 
due to the very soft form factor resulting from these potentials. Therefore we want to look 
now at the 7rN form factors extracted from the models applied in the present study. The 
values for i^jvAr(O) are also given in Table |. Note that model 2' yields the hardest 7rNN 
form factor which corresponds to a monopole form factor with a cutoff mass of about 708 
MeV whereas the model PJ of SA corresponds to a monopole form factor with a cutoff 
mass of just 317 MeV. From comparing different columns of Table | one might conclude 
that there is some influence of the form factor on the magnitude of the TBF. However, the 
variations in the individual partial wave contributions are only around a factor of two or 
three and definitely not two orders of magnitude. Furthermore, the cancellation effects are 
independent of the softness of the ttNN form factor and they tend to reduce the variations 
in the total contribution. In fact, the results for the ttN models considered in the present 
paper lie all within a range of 15 keV, as can be seen from Table |. 

At this point one may wonder how reliable results based on a rank 1 separable approxi- 
mation of the Julich models are. In order to estimate the uncertainty due to the simplicity of 
the representation we constructed another rank 1 represention where the expansion energy 
was chosen at ttN threshold (for the S-waves) or slightly above. Specifically we chose 1077 
MeV for Sn, 1000 MeV for S31 and 1100 MeV for the other partial waves. The results 
for these alternative separable representations are compared with the ones obtained for our 
'standard' choice in Table [TT]. From this Table we see that there are variations of the order 
of 20 ~ 30% in some partial waves - but qualitatively there is no change in the results. 
Therefore we are confident that the used rank-1 separable representations are sufficiently 
accurate for the aim of the present investigation. 

V. DISCUSSION 

We conclude from the previous section that the softness of the nNN form factor is not 
responsible for the smallness of the 7T7r-exchange TBF based on ttN potential models. If so, 
what makes the contribution to the three-body binding energy so small? In order to shed 
some light on this let us examine the basic two differences between the TM TBF and the one 
derived from a ttN potential model. The first difference concerns the energy dependence. 
The original ttN amplitude on which the TM TBF is based, is given in a covariant form and 
therefore depends on the energy (of the pion). However, in order to make this TBF suitable 
for application in standard (non-relativistic) 3N calculations the ttN amplitude is expanded 
in powers of k n /rriN and only the lowest order terms are kept. As a consequene the initial 
energy dependence drops completely out of the resulting TBF. In the approach of SA the 
TBF is energy dependent in a two-fold way, namely via the ttN amplitude but also via the 
nNN form factor. The effect of switching off this energy dependence in the TBF has been 
analyzed thoroughly in Ref. |l]|| where it was found that it leads to a sizeable increase in the 
resulting binding energy. However, it was concluded by SA that the approximation of fixing 
the energy in the ttN amplitude and the tcNN form factor does not lead to a sufficiently 
large change in the contributions to make them comparable with the result for the TM 
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TBF. None the less we would like to look at this point again because (unlike the potentials 
employed by SA) now the ttN interaction model itself is energy dependent and therefore the 
effects from the energy dependence might be stronger. Note that in the present case the 
energy dependence of the TBF enters at three levels: (a) the energy dependence of the ttN 
potential itself; (b) the energy dependence of the ttN t-matrix; (c) the energy dependence 
of the irNN vertex function. 

We investigate the role of the energy dependence by fixing the energy at E — 930 MeV. 
This is the value chosen for constructing the separable interaction via the EST method 
and, accordingly, where the ttN amplitudes generated by the Jiilich model and its separable 
representation agree almost exactly. Thus, the results for our rank-1 separable interaction 
should be practically identical to the one for the original Jiilich model in the particular 
case where the whole energy dependence (a)-(c) is fixed. Corresponding values are given in 
the last column of Table |T|. The numbers in the third column of Table |T| are obtained 



by fixing only the energy dependence of the nN interaction. For the ease of comparison 
results without restricting the energy dependence are also shown in the Table. We see that 
the contribution of each ttN partial wave is enhanced by about 50% after fixing the energy 
dependence of the potential. Fixing the energy dependence also in the irN t-matrix and 
the ttNN form factor leads to a further increase in the contributions from the Pu and 
P33 waves and to a slight suppression in the S31. However, those approximations definitely 
do not provide any really substantial enhancement of the resulting triton binding - in line 



with the findings of SA. All results shown in Table |T| are obtained by using model 2'. The 
other models behave qualitatively very similar and therefore we don't give the corresponding 
numbers here. 

The other major difference between the TM force and the TBF based on a ttN potential 
concerns the off-shell extrapolation of the nN amplitude. In the Jiilich model the off-shell 
properties of the ttN amplitude are completely determined by the dynamical ingredients of 
the nN model and the fit to the nN data. Since also for the potential models employed by 
SA the off-shell properties are, in principle, constrained by a fit to nN data let us emphasize 
the main difference here. In case of phenomenological separable interactions the off-shell 
behavior is, to a large extend, determined by the specific choice of the form-factor function. 
Moreover, separable interactions act only in single partial waves and accordingly the free 
parameters are determined by fitting only a single partial wave. On the other hand, in a 
meson-exchange potential like the Jiilich model the dynamical ingredients give, in general, 
contributions to all partial waves and therefore the free parameters in this model - which 



are essentially the cutoff masses in the (baryon-baryon- meson) vertex form factors [12| - 
are much better constrained by a fit to the ttN data. Clearly also here differences in the 
dynamics and/or differences in the parametrization of the vertex form factors will lead to 
variations in the off-shell properties of the resulting nN amplitude. Indeed such differences 
exist between the models 1 and 2 (or 2') considered here. But, as we have already seen in 
the last section, they to not lead to any significant variations in the results for the TBF. 
The nN amplitude used in the TM force is given by (cf., e.g., Ref. |HJ2*5|) 

T l J N {K, K) = F nNN (K)F„N N (^) {ft [a + b K ■ K + c (kl + k; 2 )] - d e^ k r k a ■ K x k^} 

(5.1) 

where i,j are pion (cartesian) indices, k^ and k^. are the momenta of the incoming and 
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outgoing (off-shell) pions, and a, b, c, and d are constants denned, e.g., in Ref. p5| . Evidently 



the off-shell extrapolation is provided by the form factor function F^nn [24] . In the standard 



version of the TM force this form factors are assumed to be of monopole type, 

F ^nn(K) = (5-2) 
K + K 

with a cutoff mass = 5.8m n ~ 800 MeV. We would like to emphasize at this point that, 
in principle, there is no connection between the ttNN vertex (with an off-shell pion) and 
the it N amplitude entering into the TBF. It would appear only in the contribution of the 
direct nucleon pole diagram to the ttN amplitude, which, however, is omitted in order to 
avoid double counting (cf. the discussion in section II). Therefore the prescription for the 
off-shell extrapolation used in the TM force must be considered as rather arbitrary. 

In the discussion above we have tacitly ignored a conceptional subtlety when we talk 
about "off-shell". In case of the Jiilich model the ttN amplitude is obtained off-energy-shell, 
while an off-mass-shell ttN amplitude is used in the TM potential. This means, that we can 
not make a simple comparison between them. The off-pion-mass-shell ttN amplitude of the 
TM force depends on the pion momenta k„- and k^.. The off-energy-shell value of the ttN 
amplitude for the potential model (at a certain energy) is given as a function of the relative 
momentum between the pion and the nucleon, Q. Within non-relativistic kinematics Q is 
obtained by (cf. Eq. flOj)) 

K - ^k N 

mjv 

where k]\r is the nucleon momentum. One can see from this relation that Q becomes equiv- 
alent to h n only in the limit of m n /m^ — > 0. 

In the following we want to discuss the off-shell properties entering into the calculations 
with the TM force and into the results presented in this paper. In view of the aforementioned 
difficulties it should be clear that any comparision can be only of qualitative nature. None 
the less, as we will see below such an analysis is useful and we believe that it indicates the 
source of the large discrepancy found in the contributions to the binding energy from the 
two approaches. 

For this purpose let us introduce a half-off-shell function in the following way, 

t in\ T a (Q, Q'; Z) 

fa(Q) = (5.4) 

where T a (Q, Q'; Z) is the off-shell ttN t-matrix (projected on the partial wave a with angular 
momentum I) at a fixed energy Z and a fixed momentum Q 1 . The factor Q l is taken out for 
convenience because then we can normalize these half-off- shell functions to 1 at Q = for 
s- as well p-waves and we can easily compare them with each other. Corresponding results 
for the Jiilich tcN potential 2' are shown in Fig. |[ where Z and Q' in Eq. ( |5.4| ) have been 
fixed to 930 MeV and 130 MeV/c, respectively. 

Let us first take a look at the p-waves and in particular at the Pu and P33 partial waves 
which provide the main attractive contributions to the TBF (cf. Table |). In this case the 
momentum dependence of the TM ttN amplitudes is roughly given by a monopole type 
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function F(Q 2 ) = A 2 /(A 2 + Q 2 ) with A = 800 MeV, cf. Eqs. (Q-jjD, which is shown 



by the dashed curve in Fig. H We observe that the corresponding half-off-shell functions 
of the Jiilich irN model fall off much faster with increasing (off-shell) momentum than this 
function. Accordingly we expect that a TBF based on the potential model will yield a much 
smaller attractive contribution to the three-body binding than one with off-shell properties 
similar to the monopole type function. 

In case of the s-waves we get large repulsive contributions from the S31 partial waves and 
small attractive contributions from S n (cf. Table |I[). Please note that also the corresponding 
half-off-shell functions are radically different. The one for S31 exhibits a strong enhancement 
whereas the one for S n falls off very strongly with increasing (off-shell) momentum. The 
s-wave part of the TM force (a- and c-terms), on the other hand, has no isospin dependence, 
cf. Eq. ( |5.1| ). This means that here the S31 and Sn partial waves have exactly the same 
momentum dependence. Therefore we suspect that a strong cancellation between the con- 
tributions from those two s-waves takes place. Indeed, actual triton calculations employing 
the TM force confirm that the total s-wave contributions are comparably small f26|| . Ev- 



idently, such a cancellation does not occur with the TBF based on the Jiilich ttN model 
because of the differences in the off-shell properties. As a consequence, the (large) repulsive 
contribution of the S31 partial wave to the three-body binding survives (cf. Table [I]). 

Summarizing this phenomenological discussion of the off-shell properties we expect that a 
calculation based on the off-shell extrapolation used in the TM force should lead to a strong 
enhancement of the attractive contributions and at the same time reduce the repulsive 
contributions. We believe that this is the basic mechanism which makes the binding energy 
obtained with the TM TBF so large. We would like to substantiate this claim quantitatively 
with a model calculation. We can do this by substituting the off-shell properties of our 
ttN model by the ones used in the TM force. This can be easily done for the separable 
representations that we are using. We only need to replace the Jiilich off-shell behavior by 
defining the form factor of the separable potential, g a (Q) '■— (Q\ V a \ij)E^), as follows: 



9a{Q) ^^a^)TTQW- (5 ' 5) 

Furthermore we fix the energy dependence in the irN amplitude and the ttNN vertex func- 
tion again. 

We demonstrate the effect on the binding energy for several values of A in Table [TV| . 
From this Table it is clear that we can get a substantial increase in the triton binding energy 
by choosing A « 600 ~ 800 MeV. Specifically one can see that the repulsion of the S31 
partial wave is suppressed and, moreover, cancels to a large extent with the Sn. At the 
same time the attraction provided by the Pn and P33 partial waves is strongly enhanced 
- as we expected from analyzing Fig. [8[ In fact, if we take into account that our results 
are based on first-order perturbation theory and therefore may underestimate the correct 
values by a factor two or even three then our simulation with the choice A = 800 MeV 



practically reproduces the TM result, which is likewise based on A w = 800 MeV. 
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VI. SUMMARY 



In the present paper we have re-investigated the origin of the large discrepancy in the 
contribution of a 7T7r-exchange TBF found by Saito and Afnan to the commenly accepted 
values obtained with the Tucson-Melbourne or Brazil TBF. Unlike SA, who employed phe- 
nomenological separable potentials, we started out from a meson-theoretical irN interaction 
model developed recently by the Julich group. This model provides a good description of 
elastic ttN scattering data. It is also in agreement with empirical information on the ttN 
amplitude in the subthreshold region. In particular, it predicts the ttN E term close to 
the empirical value. Furthermore, the decrease in the ttNN form factors F^^^{q 2 ) from 
q 2 = m 2 to q 2 = of about 4-7% and is comparable to that of the TM potential. Thus the 
form factors are much harder than those used by SA (which show a decrease of up to 20 %). 
Accordingly, the Julich irN model does not show the deficiencies which, so far, have been 
thought to be the main reason for the small contribution of the 7r7r-exchange TBF in SA's 
work. 

None the less it turned out that also the 7T7r-exchange TBF based on the 7rN amplitude 
of the Julich model is very small. The contributions to the triton binding energy are in the 
order of a few keV, which means comparable to the results obtained by Saito and Afnan. 

A detailed analysis of the main differences between the TM TBF and that derived from 
the Julich nN potential model suggests that the differences in the contribution of the TBF to 
the 3N binding energy is due to the off-shell behaviour of the non-pole ttN amplitude. In the 
Julich model the off-shell properties of the ttN amplitude are determined by the dynamical 
ingredients of the model and the fit to the nN data. As a result, the off-shell properties 
of the amplitude are different in the individual ttN partial waves. In particular, the ttN 
amplitudes that provide attractive contributions to the three-body binding (Pq, P33) fall 
off relatively fast with increasing off-shell momentum while the repulsive S31 partial wave is 
enhanced. As a consequence, the total contribution of the TBF is very small as a result of 
the cancellation effects. On the other hand, in the irN amplitude underlying the TM TBF, 
the off-shell extrapolation is done in terms of a monopole form factor with a cut-off mass of 
800 MeV for all partial waves. Since this monopole form factor falls off much slower then 
the Pu and P 33 amplitudes of the Julich nN model, the corresponding attraction provided 
by the TM force is considerably enhanced. At the same time, the repulsion in the S31 partial 
wave (and therefore any cancellation effects) is strongly suppressed. These combined effects 
do indeed explain the two-orders-of-magnitude discrepancy in the resulting contribution to 
the triton binding energy as we have demonstrated in a numerical model study. 

Naturally the question arises how realistic and well-defined the off-shell properties of the 
Julich irN model are, especially in view of the so-called quasi-potential ambiguity |27|. This 
is a topic which needs to be further investigated in the future, but it is certainly beyond 
the scope of the present paper. Here we only want to point to the fact that rather different 
ansatzes for the irN interaction (meson-exchange and simple separable forms, respectively) 
lead to qualitatively similar off-shell features and, in consequence, to similar results for the 
TBF. Qualitatively similar off-shell properties seem to be also predicted by other nN models 
- at least as far as we can judge from corresponding publications |28| , |29l . 

With regard to the off-shell extrapolation used in the irN amplitude on which the TM 
TBF is based the situation is, in our opinion, much less clear. First of all, we do not see 
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any stringent physical reason for adopting the irNN form factor for this purpose, especially 
because no 7iNN vertex is present at all in the non-nucleon-pole part of the nN amplitude 
that enters into the derivation of the TBF. Furthermore, the choice of having the same 
off-shell properties in all tcN partial waves is also hard to justify. Therefore the large 
contribution of the TM (and also the Brazil) TBF to the triton binding energy of around 1 
MeV - though certainly desired by phenomenology - must be interpreted with caution. 
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Appendix: Separable expansion of a potential with two terms 



We consider to represent a potential V which consists of two terms, V = V\ + V 2} in 
separable form. At first we expand V 2 by means of the standard EST-method, cf. Eqs. ( |2.2|) 
to Q2.6Q . Then the t-matrix T 2 (E) obtained from the separable representation V 2 agrees (on- 
as well as half-off- shell) with T 2 (E) corresponding to V 2 at the choosen expansion energies 
E = E U i = l,...,N. 

Next we expand the potential V by assuming the following form 

V(E) = \v )\ 1 (E)(v \+V 2 . (6.1) 
We determine the "form factor" \vq) in such a way that V satisfies 

V\iJ> e )=V(e)\1> e ), (6.2) 
where \ip £ ) is a solution of the scattering equation, 

\i/> e ) = \k e ) + G (e)V\il; e ), (6.3) 

at a fixed predetermined energy e. We want to emphasize that we may choose the energy e 
at which the potential V is expanded to be different from any of the energies Ei chosen for 
the separable expansion of V 2 . We see easily that Eq. ( |6.2| ) is satisfied if we choose \vq) to 
be 



and Ai(e) to be 



\v ) = {V-V 2 )\4, e ), (6.4) 



XAe) = = (6.5) 

mv-v 2 \^ £ ) 



(6.6) 



Note that the (half-off-shell) t-matrix T(E) obtained from V(e) is identical to T(E) 
obtained from V at E = e because of condition (|6.2j). 

If T(E) has a pole at E p we may choose the expansion energy e for the separable repre- 
sentation to be e = E p . The wave function at the pole, \ip p ), is a solution of the equation 

\1> P ) = G (E P )V\1> P ) . (6.7) 

The t-matrix T(E) for the separable potential V(E) is given by 

where 

E(E) = (v \G (E)\v(E)), (6.9) 
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and 

\v(E)) = (l + f 2 (E)G (E))\v ) . (6.10) 
The "form factor" \vq) is defined by 

\v Q ) = {V-V 2 )\il> p ). (6.11) 
Substituting Eqs. (|6.10|) and (|6.11j) into Eq. ( |6.9|) leads to 

t(E) = (t> |(l + G (E)f 2 (E))G (E)V - G (E)f 2 (E)\ij p ) . (6.12) 
We can evaluate ti(E) at E = E v by using Eq. ( |6.7| ), 



t{E p ) = (f \i/> p ) 
1 



\i(E p 



(6.13) 



where we have utilized, in addition, Eq. (|6.6|) . By substituting this result into Eq. ( |6.8| ). 
we see that the t-matrix T(E) obtained for the separable potential V has the same pole 
position as the t-matrix of the original potentil V. 

Let us now assume that the first term of the original potential V has a separable form, 

Vx(E) = \v )A(E)(v \, (6.14) 

where A(E) is 

A(E) = - J— . (6.15) 
E — mo 

This is exactly the case for the (direct) nucleon pole contribution, Fig. |2|(a). Then \v ) 
corresponds to the bare ttNN vertex function and uiq is the bare nucleon mass. The solution 
of the Lippmann-Schwinger equation for V = V\ + V 2 can then be written as 

T{E) = ^ E) \/ME)-UE) {V{E » + T2(£) ' (6 - 16) 

where T 2 (E) is a solution of the scattering equation for the potential V 2 . The self-energy 
£(.E) is given by 

E(E) = (v \G (E)\v(E)), (6.17) 
and the dressed ttNN vertex function, \v(E)), by 

\v(E)) = (l + T 2 (E)G (E))\v ). (6.18) 
Assuming that the t-matrix of Eq. ( |6.16| ) has a pole at E = mjv, S(mAr) is evaluated as 

Y,(m N ) 



A(m N ) 

m N — m • (6.19) 
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Also, the pole part of the t-matrix can be written in the form, 

T ^ = = l^d^^l ' (6 ' 20) 

which defines the renormalized ttNN vertex function \v R (E)). At the pole \v R (E)) is given 
by (cf., e.g., Ref. §) 

lv {mN)) = (i-EK))'/ 2 ' (6 - 21) 

where 

V(m N ) = -(v{m N )\G 2 (m N )\v(m N )) . (6.22) 

In the last step we have assumed that V 2 does not dependent on the energy E. If it does 
(like in our case) then £' is determined by 

S'W = ^(v \G (E)\v(E)) \ E=mN . (6.23) 

Furthermore, we note that the renormalized vertex function is related to the ttNN coupling 
constant f n NN and the 71NN form factor F nNN by 

W^7Viv(q 2 )a(g) = (q\v R (m N )) . (6.24) 

where a(q) is a kinematical factor depending on the particular (ps or pv) tcNN coupling (cf. 
Eqs. (3) and (16) of Ref. (|TJ). 

The wave function at the pole energy is obtained by solving Eq. 



\4> P ) = G (m N ) \v(m N )). (6.25) 

By substituing Eq. ( p.25| ) into the Ihs of Eq. (|6.7| ), we get another relation between \v{m^)) 
and \ip p ), 

\v{m N )) = V(m N )\i/) p ). (6.26) 

Now we expand the potential model V(E) using the method described above. Then the 
resulting t-matrix obtained from the separable representation has the same pole position as 
the original interaction model. We will examine whether also the 7rNN coupling constant 
and form factor determined from the separable representation agree with the ones of the 
original interaction. 

At first, we derive a relation between |-y(rajv)) and \v(tiin)). Using Eqs. (|6.10|) , ( |6.11|) , 
and finally (|6.26 ) we obtain 

\v(m N )) = (1 + f 2 {m N )Go(m N ))(V (m N ) - V 2 )\^ p ) 
= V(m N )\ip p ) 

= \v(m N )). (6.27) 
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Eq. (|6.27|) implies that the momentum dependence of the nNN vertex function obtained 
from the separable representation is identical to the one of the original potential. 

The renormalized nNN vertex function at the pole, \v R (m^)), for the separable repre- 
sentation V can be calculated in the same way as for the original potential. It is obtained 
by placing tildes over all quantities in Eq. ( |6.21| ) and ( |6.22j ): 

\v R (m N )) = " (mjv)) - , (6.28) 
' V NJ/ (1 - S'(miv)) 1/2 

t\m N ) = -(v(m N )\G 2 (m N )\v(m N )). (6.29) 
Because of Eq. |T27|) we get 

S'(mjv) = S'(mAr) . (6.30) 
Finally, by substituting Eqs. (1071) and (I6T30D into Eq. (|6T28|) it follows that 



\v H (m N )} = \v K (m N )} , (6.31) 

and consequently that f n NN = fnNN and F^nn^ 2 ) = F^nn^ 2 ). 

In the above discussion we have not specified the form of Xi(E), since we needed only 
the value of \\(E) at the energy E = to get the correct half-off shell t-matrix. If we 
assume \\{E) to be of the form 1/X^E) = E — rh , then m can be determined by Eq. (|6.6|) 
(or Eq. ( EH3D ) to be 

m = m N - (ip p \V - V 2 \ip p ). (6.32) 

Therefore, by following the outlined procedure it is possible to construct a separable 
representation where (i) the half-off- shell behavior of the non-pole part of the t-matrix (at 
selected energies), (ii) the pole position, (iii) the nNN form factor, and (iv) the value of the 
nNN coupling constant are the same as in the original interaction model. 
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TABLES 



TABLE I. The contribution of the 7T7r-exchange three-body force to the triton binding energy 



in keV for the ttN models 1, 2 and 2' of the Jiilich group HQ and PJ of Ref - @- The nNN 
form factor F w nn(q 2 = 0) determined consistenly from these models is given in the second line. 







model 1' 


model 2 


model 2' 


PJ 




0.934 


0.925 


0.962 


0.811 


Sn 


-3.8 


-2.8 


-2.3 


-4.8 


S31 


35.1 


64.4 


51.7 


26.4 


Pll 


-11.8 


-26.4 


-21.3 


-8.8 


P31 


1.8 


9.0 


6.4 


-3.6 


Pl3 


-2.2 


-0.9 


-1.1 


4.5 




-17.1 


-31.6 


-24.2 


-16.0 


Total 


2.0 


11.7 


9.2 


-2.3 



TABLE II. The effect of choosing different expansion energies for the rank-1 separable rep- 
resentation in the EST expansion. Columns labelled with Ei < E t h correspond to our standard 
choice of E$ = 930 MeV. The results given in the columns labelled with Ei > E t ^ are obtained for 
the value of 1077 MeV for Sn, 1000 MeV for S 3 i and 1100 MeV for the other partial waves. All 
binding energies are given in keV. 





model 1 

Ei > E t h 


Ei < Eth 


model 2' 

Ei > Eth 


Ei < Eth 


Sn 


-5.4 


-3.8 


-5.5 


-2.3 


S31 


37.5 


35.1 


59.0 


51.7 


Pn 


-11.8 


-11.8 


-20.6 


-21.3 


P31 


-1.0 


1.8 


0.3 


6.4 


Pl3 


-2.1 


-2.2 


-1.5 


-1.1 


P33 


-14.7 


-17.1 


-20.5 


-24.2 


Total 


2.5 


2.0 


11.2 


9.2 
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TABLE III. The effect of removing the energy dependence in the irN potential (i), and in 
addition in the irNN vertex function and the irN ampliude (ii). All results are obtained by using 
model 2'. The binding energies are in keV. 





exact 


(i) 


(ii) 


Su 


-2.3 


-3.4 


-3.9 


S31 


51.7 


76.5 


74.7 


Pll 


-21.3 


-33.1 


-39.9 


P31 


6.4 


13.7 


14.4 


Pl3 


-1.1 


0.9 


1.1 


P33 


-24.2 


-41.9 


-45.7 



Total 9.2 12.7 -0.7 



TABLE IV. The effect of replacing the off-shell properties of our ttN amplitude by monopole 
type functions (cf. Eq. |5.5| ) with different cutoff masses A. Starting point is model 2'. The binding 
energies are in keV. 



A{MeV) 


200 


400 


600 


800 


Sn 


-0.7 


-4.8 


-8.6 


-11.0 


S31 


0.4 


7.0 


13.4 


17.5 


P11 


-7.1 


-28.5 


-63.6 


-100.2 


P31 


-2.1 


3.2 


12.4 


19.6 


Pl3 


-1.1 


1.1 


5.6 


9.3 


P33 


-0.5 


-48.7 


-129.7 


-202.7 


Total 


-11.1 


-70.7 


-170.5 


-267.5 
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FIGURES 



FIG. 1. The contribution to the three-nucleon force. 



FIG. 2. Diagrams included in the ttN potential. 



FIG. 3. Comparison of the background isoscalar ttN amplitude F + (0,t). The solid 
(long-dashed) lines are obtained from the Jiilich ttN model 1 (2) of Ref. jl2j . The dash-dotted 
(short-dashed) lines are the predictions of the amplitudes that form the basis of the Tuc- 
son-Melbourne (Brazil) 7T7r-exchange three-body force. Empirical values at the Cheng-Dashen 
(CD) point, t = m 2 - and t = are taken from Ref. Q. 



FIG. 4. irNN vertex functions as function of the pion momentum in the 7riV cm. system. 
The solid (long-dashed) line denotes the prediction resulting from model 1 (2) of Ref. [12], the 
dash-dotted (short-dashed) line denotes the prediction of model 1' (2') of Ref. KM, the dotted line 
shows the result of model PJ of Ref. |0] . 



FIG. 5. Off-shell ttN transition amplitude T a (q,q'; Z) for q = q' = 130 MeV/c as a function 



of the total energy Z. The full line is the result of model 2' of Ref. 14] whereas the dashed line 



corresponds to the rank-1 separable representation described in the text. 

FIG. 6. Off-shell ttN transition amplitude T a (q,q'; Z) for q' = 130 MeV/c at three different 
energies Z as a function of the off-shell momentum q. Same description as in Fig. ||. 

FIG. 7. Illustration of the three-body force. See the text for the definitions of functions, 
momenta and energies. 



FIG. 8. Comparison of the monopole type function F(q) = A 2 /(A 2 + q 2 ) with A = 800 MeV 
with the off-energy shell behaviour of our model 2' for various ttN partial waves. f a (q) is defined 
in Eq. (Kty . 
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FIG. 2. (T.-Y. Saito et al.) 




FIG. 3. (T.-Y. Saito et al.) 




FIG. 4. (T.-Y. Saito et al.) 
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FIG. 6. (T.-Y. Saito et al.) 
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